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Rossen I. Ivanov0 
Abstract 

The squared eigenfunctions of the spectral problem associated with the 
CamassaHolm (CH) equation represent a complete basis of functions, 
which helps to describe the inverse scattering transform for the CH hi- 
erarchy as a generalized Fourier transform (GFT). All the fundamental 
0^ ■ properties of the CH equation, such as the integrals of motion, the de- 

I scription of the equations of the whole hierarchy, and their Hamiltonian 

structures, can be naturally expressed using the completeness relation 
and the recursion operator, whose eigenfunctions are the squared so- 
lutions. Using the GFT, we explicitly describe some members of the 
CH hierarchy, including integrable deformations for the CH equation. 
We also show that solutions of some (1 + 2) - dimensional members 
I of the CH hierarchy can be constructed using results for the inverse 

scattering transform for the CH equation. We give an example of the 
peakon solution of one such equation. 
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1 Introduction 

The Camassa-Holm (CH) equation [1] became famous as a model in the 
^SJ . theory of water waves. It is also known that it describes axially symmetric 

(N ■ waves in a hyperelastic rod [21 [3]. The most prominent representative of 

QO ■ the water-wave equations, the Korteweg-de Vries (KdV) equation does not 

describe the wave-braking phenomenon. In addition to the stable soliton so- 
I lutions, the CH equation, together with another recently derived nonlinear 

integrable equation, the Degasperis-Procesi equation, has smooth solutions 
that develop singularities in finite time via a process that captures the fea- 
tures of the breaking waves: the solution remains bounded, but the slope 
' becomes unbounded [4]. For the physical relevance of these two equations 

as models for the propagation of shallow water waves over a flat bottom 
one can consult, for example [3 [H [3 [HI [9l [THl [II]- More about the physical 
applications, modifications and the type of solutions of the CH equation can 
also be found in[31ll2l[l3l[Ill[l5l[Ill[I71[l8l[l9l[20]. 
The CH equation has the form 

Ut — Uxxt + 2uJUx + 3uUx — 2UxUxx — UUxxx = 0, (1) 
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where a; is a real constant. This equation is completely integrable and admits 
a Lax pair [1] 

^xx = (^ + A(m + u;))^, (2) 



= (^-nj^, + ^^ + 7*, (3) 

where 7 is an arbitrary constant and m = u — Uxx- The CH solitary waves 
are stable solitons ifa;>0[6l[l2l[l3l[2l]or peakons if a; = [H [22l |23] . 

The KdV and CH equations can also be interpreted as geodesic flow 
equations for the respective and metrics on the Bott-Virasoro group 

[MlESlESlETlEElEi]. 

The CH equation is a bi-Hamiltonian equation, i.e. it admits two com- 
patible hamiltonian structures Ji = {2u>d + md + dm), J2 = d — [Tl [50]: 

^ 6H2[m] ^ 6Hi[m] 

mt = -J2 / = - Jl / , (4) 
dm dm 

Hi = - J mudx, (5) 

= l/(„' + „.l + 2.„^)d.. (6) 

The infinite sequence of conservation laws (multi-Hamiltonian structure) 
Hn[m\, n = 0, ±1, ±2, . . ., satisfying 

J 5Hn[m] _ J 5Hn-i[m\ 
5m 5m 
can be computed explicitly [U EU [Ml El]- 



2 Generalized Fourier transform 

The so-called recursion operator plays an important role in describing an 
integrable hierarchy . The recursion operator for the CH hierarchy is L = 
Jg^^ Jl. The eigenfunctions of the recursion operator are the squared eigen- 
functions of the CH spectral problem. For simplicity we consider the con- 
crete case where m is a Schwartz class function, w > and m(x,0) > 0. 
Then m{x,t) > for all t, e.g. see [351 [15]. It is convenient to introduce 
the notation: q = m + uj. Let fc^ = — | — Xuj, i.e. 

A basis in the space of solutions of ([2|) can be introduced: f~^{x,k) and 
/"•"(x, k). For all real A; 7^ it is fixed by its asymptotic when 2; ^ 00 [35] . 
(also see [la [361 [37]): 

lim e-^'=^/+(x,A;) = 1, (9) 
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We can introduce another basis, / (x, A;) and/ (x, k) fixed by its asymptotic 
when X — > — oo for all real k ^ 0: 



lim e'''T{x,k) = 1, (10) 



Because m(x) and u are real we find that if f~^{x, k) and / (x, A;) are solu- 
tions of ([2]) then 

f+{x~k) = f+{x,-k), and r{x,k) = r{x,-k), (11) 

are also solutions of ([2]). The squared solutions are 

F^{x,k) = {f^{x,k)f, F±(x) = F(x,iK„), (12) 

where F^{x) are related to the discrete spectrum k = iKn, where < ki < 
...< Kn< 1/2. 

Using the asymptotics ([9]) , (jlOp and the Lax equation ([2]) one can show that 



L±F±(x, A;) = (13) 
A 



where 



L± = {d^ - 1)"^ Uq{x) -2 j dxm'(x)j (14) 

J ±oo 

is the recursion operator. The inverse of this operator is also well defined. 

We introduce the notation = dx. The squared solutions ()12p 
form a complete basis in the space of the Schwartz class functions m(x), and 
y, t, can be treated as some additional parameters. Also, the Generalised 
Fourier Transform (GFT) for q and its variation over this basis is 



to 



1 = ±-L r ^^^F±(x, k)dk + f; ^RtF^ix), (15) 
q{x) 2m J_^ ujX{k) ^ u;A„ 



1 

2TTi 
N 

n=l 



ojX{k) 



F'^{x,k)dk 



LoXn 



Fn{x) + -^5k^F^{x) , (16) 



where F^{x) = ^F^(x, k)\k=iKn ■ The generalized Fourier coefficients TZ^{k), 
R^, together with the set of discrete eigenvalues, are called scattering data. 
The variation is with respect to any additional parameter, e.g. y, t. 
The equations of the CH Hierarchy can be written as 



v1 



1 



0, 



(17) 
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where Pi{z) and P2iz) are two polynomials. If i}(z) = ^|fy is a ratio of 

these two polynomials one can define 0(L±) = Pi{L±)P^'^ {L±) (provided 
P2{L±) is an invertible operator). Then ()17|) can be written in the equivalent 
form 

qt + 2qux + q^u = 0, u = ^0(L±) - l) . (18) 

Due to the completeness of the squared eigenfunctions basis, from p!7|) . 
()15p and ()16p we obtain linear differential equations for the scattering data: 

nf T ikn{x~^)n^{k) = o, (19) 

i?±,±K„f](A-i)i?± = 0, (20) 

Xn,t = 0. (21) 

The GFT for other integrable systems is derived e.g. in pSj [Ml HQI HH 



Example: We now consider the case Q{z) = a_iz ^ + + aiz (where 
aj are constants). The ()17p equation can then be rewritten as 

L±^-^^^ + (a-i + aoL± + a,Ll) (yf " l) = 0, (22) 

Taking the identities 

^-1 



L±n = -2(1 - a2)"^5i^(ug^ + 2qux 
into account, we obtain an integrable equation 



qt + ai{2qux + q^u) - —q^ —{d - ^ = 0, (23) 

which becomes the Camassa-Holm equation ([T]) with the choice oi = 1, 
ao = fl-i = 0. Therefore, (|23p can be considered as an integrable 'deformed' 
version of the CH equation. Another choice for the constants, a_i = 1, 
00 = 01 = 0, leads to the extended Dym equation [H [311 [HI- If oi = 1> 
a_i = but ao 7^ the equation is usually called Dullin-Gottwald-Holm 
Equation [H [H [13 [H] . 

The Hamiltonian of (I18p with respect to the Poisson bracket related to 
the Hamiltonian operator Ji, 

lA.B} = -J(. + .a){-8---8-)a.. (24) 
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is (see [HS]) 

Jo TTuj\{ky V ^ ^^^J K 

In ()25p the Hamiltonian is given in terms of the Scattering data. In 
general, it is not straightforward to find the corresponding expressions in 
terms of the field variable q{x) (or m{x)). For example, the Hamiltonian of 
([2311 with respect to is 

= ai/7f ^ - ^/o - ^H^f, 
11 2 " 4 

where H^^ = ^ J.!^ mudx is the first CH Hamiltonian 1^ , 

is the (-l)-st Hamiltonian for the CH equation, and the integral 

/oo 
mdx = H^^ + 2uja, (27) 
-oo 

is related to the other two CH integrals [3l] 

/OO /*oo I — 

iVq-V^fdx, a= / ( /--l)dx. 
-oo J —oo ^ \ Q ' 

3 Peakon solutions of a (1+2) - dimensional equa- 
tion from the CH hierarchy 

We consider an integrable member of the CH hierarchy with two 'time' 
variables - t and y (cf. [i5] ) 

qt + KjJxy + + + + i)qx = 0, q=Ux- Uxxx + (28) 

where w, 7 and 77 are arbitrary constants. The Lax pair for (|28p is 
^xx = + A(m + w))^', 

The Lax pair represents a non-isospectral problem. Indeed, the equation 
(j28p can be written as a compatibility condition 
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where A satisfies the relaxed condition Xt — j^^Xy = 0. But we assume that 
the spectrum is t- and y-independent in what follows, and we can generalise 
the solutions obtained for the CH equation. 
We can also write (l28l) as 



iVQ)t + [(Uy + VU. + i)^]. = 0. (29) 

Then 

d^\^)t + {Uy + i^U.,+-i)^ + l3 = 0, (30) 

where (3 is an integration constant. Further, choosing (3 = — and using 
the identities 

1, (d±\V^)y\ _ 1 



we can write ( 12811 in the form 

Taking ()32p . (llSp and (I16p into account and considering variations with 
respect to y and t we obtain linear equations for the scattering data: 

nf - ^7^± ± 2ifc(7 - 1^)7^± = o, (33) 
" i^"'^ ^ " 2!;)''"^" = °' ^^^^ 

assuming A^,* = 0. For example, if 7 = = 0, then the solution is any 
function of t + 2 Ay (with appropriate decaying properties): 

7^±(y,^) =7^±(^ + 2Ay), Rtiy,t) = Rt{t + 2Xny). (35) 

We can obtain CH equation itself ([I])forx = y,-u = f7a;,7 = r/ = 0. 

We demonstrate how to write explicit peakon solutions for this equation 
= 7] = u = 0). Until now, io was strictly positive in our considerations, 
but we can take the limit u; — > [10] which produces 'peaked' solitons and 
the equations for the scattering data should also hold in this case. We assume 
that Xk = Xk{t,y), Pk = Pk{t-,y) and introduce the notations £{x) = sign(x), 
p'j.{t,y) = etc. The ansatz that produces the A^-peakon solution for the 
CH equation can be generalised as 

N 

Q{x;t,y) = ^pfc5(x - Xk). 

k=l 
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We hence obtain 

1 ^ 

U{x;t,y) = -^p,4x-Xfc)(l-e-l---'=l), 
k=l 

1 ^ 

Uy{x;t,y) = -^[pj^e(x-xfe)(l-e-l---'=l)-pfce-l---'=lxy, 
k=i 

1 ^ 

k=l 

Using this ansatz and the identity 

f{x)6'{x - xo) = f{xo)6'{x - xo) - f'{xo)6{x - xq) 

we obtain the following system of PDEs for the quantities xi^{t,y), Pk{t,y) 
from ([28]): 

1 ^ 

= ^Y.^p'ke(^^-Xk){l-e-^^'-'^>'^)-pke~^-'-->'^x',], (36) 

k=l 
1 ^ 

fe=i 

where y) = etc. The solutions of this system can be obtained from 
the A'^-peakon solution for the CH equation [22\ I23j . where the scattering 
data are now arbitrary functions of their argument (with appropriate de- 
caying properties): ii„ = i?^(t + 2A„y). For example, if = 1, then the 
system has the form 

1 / 1 / 

xi + -jPi^i = 0' Pi + -jPiPi = 0> 

with a solution Pi = — xi" = const, xi = lni?i(t + 2Aiy), cf. [221 ES]- The 

system for = 2 peakons is (we assume that xi < X2 for all y and t, i.e. 
the case for which the A^-peakon solution for the CH equation is obtained 
in [221 [23]) 

XI = ^[-P2(l-e-'"^-"^l)-pix;-p2e-l"i-"^l4], (38) 

X2 = ^[p'i(l-e-l^^-^^l)-pie-l^^-^^lx;-p24], (39) 

Pi = -^Pi[K+P2e-l"^-"^l-p2e-l"^-"^lx^], (40) 

P2 = -^P2[p'ie-l^i-^2l+p^+pie-l^i-^2lx;], (41) 
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with solutions [22l [23] 



Ai A2(Aii?i + A2-R2) ' Aii?i + A2-R2 

(Ai - X2)^RiR2 , / J, , D \ 

Xl = In 2 2 , X2 = ln(i2i + i?2), 

where Rk = Rkit + ^X^y) (this can be easily verified). We note that the 
total momentum pi + P2 = — \^x^ is conserved. Of course, in general, 
for an arbitrary N , the A^-peakon solution for the CH equation obtained 
in |22t [23] can be used because the inverse scattering method for hierarchy 
psp is the same as for the CH equation |341 [36] . The only difference is the 
time dependence of the scattering data (and/or the additional y-dependence, 
etc.). In this example, y has the meaning of a second 'time' variable. Clearly, 
the conserved quantities in terms of Xk and pk have the same form as those 
for the CH peakons and can be expressed in terms of the quantities Afc, which 
we already assumed to be independent of y and t. But the Hamiltonian 
formulation is problematic because formally ^{z) = for the peakon solution 
and the Hamiltonian with respect to ([M|) is degenerate because of ([^S]) . 
Moreover, the right-hand side of system (I38p - ()4ip involves not only the 
quantities x^ and pk but also their y-derivatives. 

The explicit dependence on the scattering data given in [36j can be used 
in the same way for the A'^-soliton solutions of the CH hierarchy. The sit- 
uation where the initial data condition g(x,0) = m(x,0) + uj > Q does not 
hold is more complicated and requires a separate analysis 01 [351 [13 • 
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